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Abstract. In real- world sampling of k-space data, one generally makes a stochastic 
error not only in the value of the sample but in the effective position of the drawn 
sample. We refer to the latter as imprecise sampling and apply this concept to the 
fourier-based acquisition of magnetic resonance data. The analysis shows that the 
effect of such imprecisely sampled data accounts for contributions to noise, blurring, 
and intensity-bias in the image. Under general circumstances, the blur and the bias 
may depend on the scanned specimen itself. We show that for gaussian distributed 
imprecision of k- vector samples the resulting intensity inhomogeneity can be explicitly 
computed. The presented mechanism of imprecise k-space sampling (IKS) provides 
a complementary explanation for the phenomenon of central brightening in high-field 
magnetic resonance imaging. In computed experiments, we demonstrate the adequacy 
of the IKS effect for explaining central brightening. Furthermore, the experiments 
show that basic properties of IKS can in principle be inferred from real MPJ data by 
the analysis on the basis of bias fields in magnitude images and information contained 
in the phase-images. 
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1. Introduction 

In this paper, a concept we would like to refer to as imprecise sampling is introduced - 
and applied to MR image acquisition experiments. We will show that imprecise sampling 
provides an alternative explanation for the mechanism of central brightening. This may 
be of interest in cases where the common interpretation in terms of dielectric resonance 
(e.g. Bomsdorf et al 1988) may fail. That may happen when high tissue conductivity 
prevents the creation of strong dielectric resonances (Collins et al 2005). The discussion 
adds to the recently revived discussion on this topic (Collins et al 2005, Wiggins et al 
2005, Tropp 2004, Yang et al 2002, Ibrahim et al 2001, Roschmann 2000, Kangarlu et 
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al 1999). Nonetheless, the concept of imprecise sampling can as well be discussed in a 
self contained fashion. 

Imprecise sampling basically relies on the observation that in real-world 
experimental situations sampling some property with respect to a continuous domain 
of control-parameters, such as sampling the k-space in an MRI experiment, can not be 
performed with infinite precision. We may intend to experimentally prepare and then 
sample the property for a specific value of the control parameter, i.e. k, the index of the 
fourier-mode, but we will actually draw the sample somewhere near its intended value 
k. This is, we will sample the property at some k + 5k, where 5k is a deviation from 
the precise k-space location k, the imprecision. Since standard Magnetic Resonance 
Imaging (MRI) is based on sampling the fourier-space of the underlying spin-relaxation 
distributions of the specimen with the k-space vectors k as control-parameter of the 
sampling process, the concept can directly be applied to this experimental situation. 
We will refer to this as imprecise k-space sampling (IKS). 

We have so far stated that imprecise control on the sampled fourier-mode inflicts 
deviations <5k of the true fourier-mode location k but we have not specified any properties 
of the inflicted deviations. Still, imprecision can be interpreted as some kind of noise 
corrupting the control-parameter. In the following, we will model the imprecision 5k as 
a stochastic process and the properties of the imprecision can be expressed in terms 
of properties of the stochastic process. These properties have to be specified and 
argued. For instance, in a multi-echo sequence the fourier-modes are sampled along 
the read-direction by subsequent refocusing pulses and switching of read-gradients. The 
acquired samples will subsequently get more imprecise with each refocusing pulse and 
the cumulative error can be seen as a Wiener process. Similar arguments may hold in 
different sampling strategies as for instance spiral acquisition procedures. However, the 
contributions of this noise to the resulting image go beyond the additive noise model 
which is well known in the context of linear degradation models (LDM). 

In many situations, image post-processing relies on the observation that there 
are several sources of degradation of an image that usually are classified into several 
categories, e.g. blur, noise, and intensity inhomogeneities (bias). These categories 
can be cast into degradation operators providing building blocks for LDMs which have 
successfully been utilized in image reconstruction and restoration methodology (Geman 
and Yang 1995, Husse and Goussard 2004). However, when nonlinear effects of k-space 
imprecision become relevant, we can show that IKS provides a unified approach for 
explaining contributions to blurring and intensity inhomogeneities. 

The theory of MR bias-fields is usually based on the analysis of the electro-magnetic 
field properties under specific geometric conditions and dielectric properties of the 
specimen in the scanner, e.g. (Hoult 2000). For known geometries, the Maxwell- 
equations can be computed and bias fields predicted, e.g. (Wang et al 2002, Jin 1996). 
Usually, the precise imaging geometry is not available and image post-processing is 
confronted with the challenge of correcting intensity-inhomogeneities for instance for 
de-biasing, segmentation and tissue classification (Ahmed et al 2002, Prastawa et al 
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2005, Luo et al 2005). All these explanations or predictions of bias-fields in general and 
central brightening in particular argue very closely to the level of physics, i.e. electro- 
magnetism. 

Central brightening is a well-known phenomenon in the context of high-field MRI 
(Collins et al 2005, Tropp 2004). When the wavelength of the RF pulse used for 
slice selection, i.e. nuclear spin excitation, gets comparable with the dimension of the 
specimen the dielectric properties of the specimen, i.e. its tissue types, may influence the 
high frequency waves. This, for example, has been pointed out for cylindrical resonators 
by Tropp (2004). Alternative explanations have been sought when dielectric resonance is 
questionable and, for instance, Collins et al (2005) argues that high tissue conductivity 
prevents strong dielectric resonance. It is interesting to note that according to collected 
measurements of mammal tissue-conductivity (e.g. Gabriel et al 1996) changing the 
frequency domain from 127 MHz (3T) to 380 MHz (9T), white- matter renders an 
increase of tissue conductivity by a factor 2-3. However, central brightening is a sub- 
topic of the significant intensity inhomogeneities observed in high-field MRI (Vaughan 
et al 2001, Bomsdorf et al 1988) posing not only challenges for optimal coil design (e.g. 
Vaughan et al 1994) but also for the precise understanding of physical and physiological 
processes under the conditions of the given context. 

With IKS, we take a complementary approach to this phenomenon respecting 
the physical level, in terms of explaining the finite control on experimental control 
parameters, but dealing with the problem on the level of data acquired under conditions 
of imprecise control. We show later that for specific reasonable properties of the k-space 
imprecision <5k the intensity-inhomogeneity from IKS can be computed explicitly and 
implies a gaussian shaped central brightening of the fourier-reconstructed images. The 
purpose of this paper is to present and demonstrate the concept of IKS. 

The paper is organized as follows. We begin by introducing the IKS concept 
with respect to fourier-encoding MRI and how the imprecision can be seen as noisy 
control-parameters k. We continue by showing how k-vector imprecision can be made 
responsible for contributions to visible blur, bias, and noise. We present an example 
that can be computed explicitly. Moreover, this example reproduces bias-fields with the 
characteristics of central brightening. We close by presenting characteristic computed 
images of IKS and discuss the results. 

2. Fourier-encoded MRI and noise 

In order to apply the idea of IKS to MRI, we recall that fourier-encoding is achieved 
by imposing phase-differences on the excited spins rotating at Larmor frequency by 
switching on linear magnetic gradients B«(x) = <?jXj, where gi are the constant gradient 
strength values, in {% = 0) phase- and (i — 1) read-direction. Without going into 
details, we can say that the phase-differences 84>i oc Bj are proportional to the magnetic 
field strength in the respective directions. The trick of fourier-encoding is to choose all 
effectively involved parameters such that the phase-shifts encode a sample in k-space, 
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i.e. J2i 34>i — kx where k is the ideal fourier- vector, x is indexing the field of view 



The k-space vectors k are indeed control parameters of an MR imaging experiment. 
It is obvious that we can not achieve control with infinite precision. First of all, 
an MR scanner is a real device with finite accuracy in triggering read-out times and 
switching phase- and frequency-encoding gradients. Secondly, tissue may have different 
properties of interacting with the magnetic field related to the dielectric constant or the 
conductivity of the media. Now, suppose that the dielectric resonance is weak due to 
high tissue conductivity. Then, switching magnetic fields will induce stronger (ionic) 
currents in the tissue which adds to the effective magnetic field and locally distorts the 
magnetic field strength controlling the relative phase-shifts of the rotating nuclear spins. 
Yet, exactly these phase-shifts are preparing the experimental fourier- vector associated 
with the sampled fourier-mode. This tells us that the overall effect of these influences 
can be cast into an imprecision of the control-parameter, i.e. k — > k + 5k, with the 
effect that we do not measure the ideal fourier-modes of the sample-distribution but 
some effective fourier-mode 



where p(x) is the ideal, static, spatial distribution of the spin-relaxation. 

When we assume conditions where the interaction of the specimen and the magnetic 
field of the MR scanner lead to stochastic variations of the local magnetic field strength 
within the specimen, then the contributions to imprecision 5k for a specific k have to 
be drawn from all over the specimen, i.e. for all locations x, so that 5k = 5k(k, x) is 
depending on both k and x. In this sense, 5k functionally depends on the spin-relaxation 
distribution p of the specimen. Different tissue classes, accounting for possibly non- 
uniform spatial characteristics, may cause different local characteristics of imprecision 
for different values of the control-parameter k. 

3. Bias, Blurring, and Noise 

In this section, we discuss how k-vector imprecisions adds to degradations, i.e. to 
intensity inhomogeneities, blur, and noise, of the fourier-reconstructed images. Let () 
denote the usual ensemble average or expectation-value with respect to the stochastic 
process 5k. Defining 



the expected effect of gradient noise on the fourier-basis vectors is 

(exp (i(k + 5k)x)) = G(k, x) exp (ikx). As a consequence, the average image of the 

simple fourier-reconstruction appears to be 



(FOV). 




(1) 



G(k,x) = (exp (iSkx)) 



(2) 




(3) 
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where p* = (p e ff) is representing the expectation on an ensemble of recorded data. 
Further, defining 



We see that g(x' — x, x) has characteristics of a point-spread function via the x' — x 
dependency with a local condition on x which refers to a bias. In case 5k functionally 
depends on p, this dependency will be inherited by g and the observable bias and blur 
will depend on the underlying nuclear spin relaxation distributions. 

Note that, <?(x', x) basically can be a complex function. However, when — 5k(— k, x) 
is equivalent to 5k(k, ctS cl stochastic process then g is a real function and no informa- 
tion should be found in the MR phase image. Moreover, when (5kf) = for all odd n, 
then G is a real function and g(— x', x)* = <?(x', x), where * denotes complex conjugation. 

We would like to note that the derivation of equation (5) is based on the average 
data p*. Simulated single images with an imprecision 5k sampled on the image grid 
size are extremely corrupted and it takes an average over many images to produce a 
reasonable signal to noise ratio (SNR). We might ask why it is possible to see the central 
brightening effect in single real recorded images after all? To see this, one should note 
that in the MRI experiment the fourier-modes are sampled one by one. Yet, each time 
we sample some specific k, i.e. p e ff(k), the integral in equation (1) is collecting the 
stochastic events of imprecision over the entire spatial domain. This is comparable to 
the ensemble average in case the spatial correlation length of the stochastic process 5k 
is small with respect to the voxel size of the image-grid. The spatial correlation length 
defines the voxel-size of a fine grid. Imprecision simulated on the fine grid can then be 
down-sampled to the coarser image grid. But down-sampling is equivalent with averag- 
ing over a large number of coarse-grid-images whith their coarse- voxel values drawn from 
the ensemble of fine-voxel values contained in the coarse-voxel. Finite but small spatial 
correlation lengths of 5k comparable to the voxel-size, however, will not be completely 
smoothed by integration over the domain of a coarse- voxel and the events of imprecision 
remain visible after fourier-transforming back to x-space. This remaining imprecision is 
basically not distinguishable from additive noise while the contributions that have been 
effectively averaged by integration appear as blur and bias. 

Example: Let us look at the stochastic process 5k(k, x) and assume the components 
5ki to be normally distributed with zero mean and co- variance -D^k, x) = (5ki5kj). 
The standard deviation is denoted a with o l o = D, where a 1 denotes the transposed 
matrix of a. Let us further assume point-wise statistical independence of the compo- 
nents 5ki, i.e. (<5k™5k™) = (5k™) (5k") for all (k, x). This can be justified since phase- 
and frequency-encoding gradients are generally not switched on at the same time and 




(4) 




(5) 
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implies that D and thus er are diagonal. Then it is possible to explicitly compute G and 
g. Starting with the definition of G in (2) we expand the exponential into a series, and 
use that the (2n)'th central moment of a gaussian process equals to (cr 2 /2) n (2n)!/n!, 
and all odd moments vanish. This implies that in the expansion 

G(K *) = E E (") <<*ox o r m (to)™ (6) 

n=0 V 1 )- m=0 V'V 

only terms with even n and even m contribute. After minor manipulations, one arrives 
at the following expression: 

G(k,x) = exp(-i| ( 7(k,x)x| 2 ) . (7) 

Suppose a, or equivalently G, is not depending explicitly on k, then the definition of g 
in (4) implies that 

4(x'-x,x)=5(x'-x)G(x) (8) 

and thus 

|dxs(x'-x,x)p(x) = G(x')p(x') . (9) 

In this case we can not only compute the bias field explicitly but also note how the bias 
field may functionally depend on p via the x-dependency of a, as different tissue types 
correspond to different spatial domains with different standard deviations of imprecision. 
When the imprecision is k-independent no blurring is introduced in the image and we 
have pure bias. When a is not depending on x as well, i.e. there is no functional 
dependency on the ground truth, then we explicitly computed a gaussian intensity bias 
with the maximum in the image center. We want to point out how this simple example 
in particular demonstrates how IKS, as a general concept, provides a complementary 
explanation for the phenomenon of central brightening. Moreover, we have to note that 
in the general case a = <r(k, x) the bias-point-spread-function g can not be written as in 
(8) with the immediate consequence that noise levels a varying with k have to introduce 
blurring or distortions in the image of p. 

4. Experiments 

One can alternatively compute the IKS effect according to equation 5 or simulate the 
IKS effect. The simulations can be done by effectively computing fields 5k(k, x) for every 
k, multiply the associated factor exp(i5kx) point-wise with the phantom-distribution 
p perform a fast-fourier-transformation (FFT), and only choose one /c-mode according 
to the sampling sequence. After all fourier-modes have been sampled in this manner 
we transform the sampled fourier-information back to x-space with an inverse FFT. 
Simulations are rather time intensive compared to computing images but demonstrate 
the contribution of IKS to noise, which is lost in the computed images. In figure (1, b,c) 
we show this for constant a. 
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(a) Phantom (b) Simulated (c) Computed 



Figure 1. Simulated and Computed IKS effect for constant a. 

Our experiments consist of computed and simulated effects of IKS on a known 
software-phantom p(x) (figure 1, a). In order to get an impression of the IKS effect we 
look at solutions of the example given in section (3) and choose some simple models for 
the standard deviation of the imprecision 

o/(x) = o- p(x) Q , 

< r JJ (k)=(7o(l + |k|)' 3 , (10) 
om(k) = <T (l + 7(k - fc off )) , 

where do is an amplitude factor, a, (3, and 7 > supply an additional parameter to each 
model. In the am model /c Q g is an offset such that for the first sample in read direction 
(i = 0) <5k has a standard deviation a and the standard deviation of the imprecision 
increases with <7 7Ak for any subsequent fourier-mode in read direction, where Ak 
is the increment between two subsequent fourier-modes in read direction. In the first 
model we can associate different levels of imprecision with different regions of intensity- 
values of p representing the different tissue types. While 07 models imprecision that 
may vary depending only on the origin of imprecision, i.e. its functional dependency on 
p, the models an and o~m are roughly modeling effects that can be associated with the 
methodology of sampling. Particularly am can be regarded as a model for the linear 
increasing standard deviation of a Wiener-process in read direction. 
We note that these choices have been made in order to visualize some basic properties 
of the IKS effect. Our choices are extremely simple models that still have the capac- 
ity to demonstrate what effects we may expect from IKS with only a few parameters. 
Moreover, since we work with simulations on simple software phantoms, we only have 
to discriminate a small number of different "tissues". A power of p is therefore sufficient 
to explore the situation in x-space. In k-space, a power in |k| can roughly be associated 
with a spiral or radial sampling protocol, while the linear model am m &y be associated 
with the cumulative error made in preparing the fourier-modes in a multi-echo sequence 
in read-direction. 




Figure 2. Model 07, with a = — 1 (a), a = —0.5 (b), a = 0.5 (c), a = 1 (d); 
The amplitude 00 has been chosen such that the smaller circular structure is 
rendered with comparable gray values. 



The experiments show the following. First of all, that the simulated effect agrees 
well with the computed effect (figure 1). The simulated image is the average of 40 images 
with simulated imprecision. The images in figure (2) are computed with the imprecision- 
model 07, which is not depending on x. We can see how different choices of a associates 
different radial decay rates to the different "tissues", which demonstrates how the effect 
of IKS functionally depends on the specimen. Stronger tissue related imprecision causes 
a stronger radial decay of intensity in the tissues gray-values. Moreover, the imaginary 
part of these images vanishes up to computing accuracy so that one may conclude, 
that it is not the x-space related functional dependency of imprecision, that causes 
non- vanishing imaginary images. 

The images in figure (3), computed with the imprecision-model an and am, 
show that k-depending imprecision introduce non vanishing imaginary parts of the 
resulting images and therefore also non- vanishing phase-images (figure 3, b,d). The 
magnitude images (figure 3, a,c) show IKS for the radial and the spin-echo model and it 
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Figure 3. (a) amplitude- and (b) phase- image of the radial model onifi = 0.1) 
(c) amplitude- and (d) phase-image of the multi-echo model (7/77(7 = 0.001) 



is obvious that the radial imprecision introduces a much weaker phase-image although 
the amplitude images show comparable central brightening. Furthermore the procedural 
asymmetry between phase-encoding and read - direction in the spin-echo model can be 
noticed clearly. 

In consequence it is, at least in principle, possible to study imprecision on the basis 
of real experiments. The problem can be separated into contributions to the IKS effect 
from x-space, regarding to central brightening decaying at different rates in different 
tissue types, and contributions from k-space, i.e. the sampling process itself, regarding 
to information contained in phase-images. One may also note the implication that under 
specific circumstances it should be possible to distinguish tissue with identical spin- 
relaxation properties Tl and T2 and thus basically have to be rendered with identical 
gray-values as long as they introduce different levels of imprecision. In this hypothetical 
situation the tissue would become distinguishable under high-field conditions due to 
different radial decay-rates of the image intensity. Imprecision may therefore possibly 
provide information that is complementary to spin-relaxation and regards to properties 
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like the dielectric constant and the tissue conductivity. 
5. Conclusions 

We have introduced the concept of imprecise k-space sampling (IKS) and demonstrated, 
that this concept provides an alternative explanation for the phenomenon of central 
brightening in high field MRI. We have demonstrated by computer experiments, on 
the basis of simple imprecision-models, that (i) such imprecissions do have the effect the 
theory predicts, (ii) functional x-distribution dependency alters tissue-decay rates, while 
it can not be made responsible for non-vanishing imaginary-images, (iii) functional k- 
distribution dependency is capable of producing non-vanishing imaginary images. These 
results show, that it is possible to collect information on the imprecision from the 
complex fourier-data of an MRI experiment by analysis of intensity inhomogeneities 
and phase images. However, that imprecise controll has consequences beyond noise, is 
a finding as basic as surprising. 
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